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PREFAC: 


^ Several  recent  Russian  oceanographic  contributions  have  been  ~ 
translated;  they  are  germane  to  the  topics  of  ocean  dynamics,  e.g., 
turbulence  and  internal  waves . The  article  by  R.  V.  Ozmidov  on  oceanic 
turbulence  is  one  of  the  clearest  discussions  of  the  topic  to  date* and 
it  also  presents  new  quantitative  results.  The  application  of  the 
theory  of  random  function  to  ocear.  currents  by  B.  N.  3elyaev  and  U.  S. 
Boldirev  suggests  avenues  of  approach  to  the  prediction  of  ocean  currents 


where  a good  statistical  knowledge  of  the  current  field  is  available. 
The  other  articles  give  further  theoretical  and  experimental  considers 
\tions  related  to  periodic  motion  in  the  sea. 


The  intrinsic  value  of  this  material  is  twofold:  it  presents 
techniques  and  data.  Its  significance  to  underwater  acoustic  systems 
and  research  is  obviously  that  sane  of  the  variability  observed  in 
acoustic  transmission  data  may  be  attributed  to  motions  of  the  medium 
if  this  variability  is  to  be  understood  r sane  knnwi  pdprp  of  the  medium 
variability  must  be  achieved  S V /X  a 1 1 /\  A 7 
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3CME  ASPECTS  OF  OCEANIC  TURBULENCE* 


Ly 

R.  V.  Ozmidov 

Institute  of  Oceanology,  Academy  of  Sciences,  USSR 


ABSTRACT 

Semi -empirical  and  statistical  methods  for  the 
investigation  of  oceanic  turbulence  are  discussed. 

The  first  method  concerns  the  difficult  problems 
of  averaging  different  characteristic  fields  in 
the  sea  and  of  determining  virtual  exchange  coef- 
ficients. These  exchange  coefficients  depend 
critically  upon  the  scale  of  phenomena  involved. 

The  second  method  has  appeared  in  oceanographic 
practice  only  in  recent  years.  Results  are  given 
for  the  calculation  and  analysis  of  the  statistical 
characteristics  of  large-scale  turbulence  in  the 
Hack  Sea.  Data  of  current  speeds  Trot  multi -day 
stations  at  the  20  and  100  meter  horizons  are  used. 

The  following  computations  are  made:  turbulent 
intensity,  tensor  components  of  correlation  moment, 
correlation  coefficient,  autocorrelation  coefficient, 
and  spectral  and  structural  functions.  This  analysis 
leads  to  conclusions  about  the  anisotropy  of  oceanic 
turbulence  and  the  ratio  of  inertial  to  turbulent 
motions  in  the  sea.  The  amplitude  of  inertial  os- 
cillations cuid  the  rate  of  dissipation  of  turbulent 
energy  are  estimated. 

If  any  sensitive  device  for  the  measurement  of  current  speeds  is 
lowered  to  3ome  depth  in  a region  of  the  sea,  described  cm  a chart  as 
a region  of  constant  current  (directed  for  example  to  the  north),  and 
if  the  device  is  held  at  the  sample  point  for  a few  minutes,  then. 


: 


^-Published  in  Journal  of  Oceanographic  Society  of  Japan.  Twentieth 
Anniversaiy  Volume,  1962;  article  received  May  22,  1962 


after  reading  it,  there  will  be  a high  probability  that  the  device  will 
shew/  that  the  currents  are  by  no  means  directed  to  the  north  and  that 
the  absolute  magnitude  is  substantially  different  frcm  that  anticipated 
A measurement  of  speeds  at  that  same  point  in  the  sea  at  some  other 
moment  of  time  may  give  a different  result  frcm  the  former  and  may 
again  not  coincide  with  that  inferred  from  the  data  on  the  chart. 

Only  by  producing  a lengthy  series  of  observations  and  averaging  them 
may  we  hope  to  acquire  a chart  which  is  close  to  what  is  displayed  on 
our  chart.  That  example  shows  that  constant  (or  slowly  and  regularly 
changing)  currents  in  reality  do  not  exist  in  the  seas  or  oceans  at  a 
given  moment  of  time,  because  they  are  the  result  of  the  application 
of  an  averaging  operation  to  actual  instanteous  currents  in  natural 
fields  of  turbulent  velocities.  Reynolds  (0.  Reynolds,  1895)  suggested 
that  the  representation  of  components  of  such  instantaneous  turbulent 
speeds  be  in  the  form  of  the  following  sums: 


Mse-tV;  vt=0+t/;  M >mw+ui‘ 


where  * , 0 and  S>  are  mean  values,  and  i/,  i/  and  tv'  are  the 

fluctuating  (turbulent)  components  of  speed. 

Ey  "mean  speed",  an  expression  of  xhe  following  form  is  ordinarily 
understood: 


where  T is  the  averaging  period.  Similar  to  such  time  averages,  i 
integration  can  be  performed  an  eny  spatial  field  averaged  over  G 
linear  dimension. 


Obviously,  as  T increases,  more  large-scale  pulsations  of  speeds 
are  included  in  the  averaging.  If  this  averaging  period  significant!; 
surpasses  the  duration  of  all  the  flew  oscillations,  then  the  average 
value  will  not  be  a function  of  time,  and  repeated  averages  will  not 
change  their  values.  Consequently,  in  that  case,  the  following  rela- 
tionships will  be  true: 
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Such  a theory  for  averaging  fields  of  speeds  is  related  to  the 
usual  theory  of  the  distribution  and  transmission  of  energy  on  various 
scales  of  motion  in  the  ocean.  It  is  believed  that  energy  from  surface 
sources  (tangential  wind  friction,  or  nanuniform  heating  of  water 
masses  by  the  sun)  causes  the  large-scale  mean  movement,  which  in  its 
turn  supplies  energy  to  smaller  structures  in  the  flow,  and  so  on,  to 
the  smallest  turbulent  vortices,  in  which  there  occurs  a dissipation 
of  the  mechanical  energy  into  thermal  motion  in  the  liquid.  That  would 
seem  to  be  the  natural  scheme,  serving  as  the  basis  for  many  theoretical 
constructions;  however,  we  do  not  have  sufficient  experimental  confirma- 
tion concerning  conditions  in  the  ocean,  and  this  theory  has  even  been 
doubted  in  recent  years.  Stommel  (H.  Stcranei,  1949)  considered  that 
the  influx  of  energy  to  oceanic  currents  from  the  tangential  friction 
of  wind  existed  not  only  for  mean  flow  but  for  the  entire  spectrum  of 
turbulent  oscillations.  Moreover,  recently  a supposition  has  been 
expressed  about  the  possibility  of  a reverse  flow  of  energy,  i.e., 
from  large  turbulent  eddies  in  the  ocean  (meanders)  to  the  mean  current. 
Confirmation  of  such  a possibility  was  recently  introduced  by  Webster 
(F.  Webster,  1961),  who  calculated  the  transport  of  kinetic  energy 
from  large  vortices  (meanders)  to  the  mean  current  of  the  Gulf  Stream. 
This  placed  in  doubt  what  is  ordinarily  accepted  in  the  theory  of  ocean 
currents  as  an  axiomatic  statement  concerning  the  analysis  of  external 
forces  acting  on  a given  water  mass,  i.e.,  that  we  ought  to  obtain  from 
the  average  aquations  of  motion  for  fluids  some  regular  (either  station- 
ary or  regular  depending  upon  the  relevant  forces)  speed  field  in  the 
observed  region  of  the  ocean.  Realistically,  if  there  is  initially  a 
chaotic  transfer  of  turbulent  oscillations  (even  though  they  would  be 
of  a large  scale),  one  must  apply  what  is  termed  the  "average"  operation 
to  the  current  and  test  for  the  irregular  variations,  which  may  not  be 
used  for  evaluating  the  average  equations  of  motion.  In  that  case, 
numerous  theoretical  statements  concerning  oceanic  currents  and  turbulent 
exchange,  which  are  used  for  the  average  equations  of  motion,  diffusion, 
and  heat  conductivity  ought  to  be  adaptable  to  the  theoretical  deter- 
mination of  the  actual  fields  of  analyzed  values  through  the  instan- 
taneous functions.  But,  for  the  condition  of  validity  for  ordinary 
schemes  of  turbulent  oceanic  motions,  i.e.,  that  there  exists  an  energy 
balance,  the  knowledge  of  either  the  speed  or  other  average  character- 
istic fields  (temperature,  salinity,  pressure,  etc.)  is  often  found 
quite  insufficient  for  solving  many  problems.  So  for  the  prediction 
of  the  trajectory  of  motion  of  a ship  in  an  associated  current  body 
(e.g.,  for  a drifting  ship),  it  is  necessary,  obviously,  to  have 
information  about  the  characteristics  of  the  instantaneous,  specific 
fields  which  move  this  body. 


5 


USL  Tech  Memo 
No.  913-27-64 


Investigations  into  the  characteristics  of  instantaneous  turbulent 
speeds  in  the  ocean  began  to  develop  in  recent  years,  and  the  direct 
results  achieved  are  still  quite  modest.  Moreover,  expressions  des- 
cribing the  field  of  turbulent  speed  are  very  complicated,  and 
difficulties  of  a fundamental  character  arise  in  their  solution.  For 
this  reason,  the  equations  are  ordinarily  averaged  in  some  mode  and 
brought  into  concurrence  with  the  form  of  the  equations  for  laminar 
mean  flow,  and  the  effect  of  turbulence  is  accounted  for  by  the 
introduction  of  sane  virtual  coefficients  of  turbulence.  These  coef- 
ficients, which  may  be  determined  from  the  equations  of  motion,  do  not 
appear  to  be  physical  constants  and  depend  upon  the  characteristics  of 
the  flow  and,  most  important,  upon  the  period  (or  scale)  of  the  average 
used  in  each  problem.  That  consideration  does  not  have  special  signif- 
icance for  those  turbulent  flows  in  which  there  are  maximum  periods  of 
turbulent  micro-pulsations,  i.e.,  whose  averages  became  stationary 
without  further  increase  of  the  averaging  period.  Oceanic  turbulence 
is  a considerably  complex  case.  Here  it  is  difficult  to  identify  the 
value  of  T over  which  all  the  turbulent  pulsations  in  the  sea  would  be 
eliminated.  It  follows  that,  in  the  sea,  the  choice  of  averaging  period 
T (or  scale  for  a spatial  average)  plays  a determining  role  in  the 
empirically  derived  coefficients  of  turbulence.  However,  in  the  initial 
studies  of  turbulence  in  the  sea,  when  coefficients  of  turbulence  were 
not  determined  in  the  basic  measurements  of  fluctuations  and  investiga- 
tions of  magnitudes,  but  by  the  method  of  observing  averaged  fields  of 
these  quantities,  questions  about  scales  of  averages  (scales  of  phen- 
omena) did  not  occupy  a significant  influence.  Meanv/hile,  just  the 
differences  in  scales  of  averages  (phenomena)  were,  to  a significant 
degree,  the  explanation  for  the  enormous  diversity  in  the  magnitudes 
of  the  coefficients  of  turbulent  exchange  (vertical  as  well  as  horizontal) 
acquired  by  various  authors  for  the  different  regions  of  the  world  ocean. 

Some  consideration  is  required  regarding  the  above  described 
method  of  constructing  and  solving  the  problem  of  the  dynamics  of 
turbulent  motions  and  the  processes  of  turbulent  exchange,  by  wiiich 
the  semi -empiri ceil  coefficients  of  turbulence  are  determined.  It  is 
to  be  noted  that  most  of  the  major  problems  in  the  dynamics  of  sea 
currents  and  turbulent  exchange  in  the  sea  are  at  present  determined 
practically  only  in  the  limits  of  such  semi -empirical  theories.  That 
circumstance  explains  the  limitations  of  most  of  the  works  on  oceanic 
turbulent  exchange  (momentum,  heat,  salt,  etc.),  which  are  devoted  to 
various  indirect  methods  of  determining  turbulent  coefficients.  It  is 
hardly  expedient  to  enumerate  these  works  because  the  majority  of  them 
suffer  by  being  largely  insufficient  due  to  the  absence  of  data  about 
the  relationship  between  the  coefficients  of  turbulence  found  from  the 
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determined  parameters  of  each  specific  problem  and  the  scale  of  average 
used  in  their  solution.  We  come  only  briefly  to  a series  of  works  in 
which  this  insufficiency  is  absent  and  which  is  devoted  to  investiga- 
tions into  the  dependence  of  exchange  coefficients  upon  scale  of  the 
average  (of  the  phenomenon) . 


These  investigations  were  principally  made  through  observations 
of  the  motion  of  discrete  particles  in  associated  oceanic  currents. 

If  data  about  the  positions  of  a series  of  such  particles  at  successive 
moments  of  time  exist,  the  coefficient  of  turbulent  exchange  (diffusion) 
may  be  determined  by  the  formula : 


where  M is  the  change  in  distance  l (scale  of  phenomenon)  between 
individual  pairs  of  particles  occurring  in  the  time  Jt  . The  bar 
indicates  an  average  determined  over  a series  of  individual  observations 


The  first  of  these  experiments  in  the  sea  were  achieved  by  Richard 
son  and  Stommel  (L.  H,  Richardson,  H.  Stomcel,  194S;  H.  Staamel,  1949) 
and  repeated  on  a broader  interval  of  scale  change  (up  to  / = 10 ^ cm) 

by  R.  V.  Ozmidov  (1957,  1959).  These  experiments  shewed  that  the 
coefficient  of  turbulence  in  the  sea  really  depends  upon  the  scale  of 
the  phenomenon,  for  that  dependence  is  expressed  as  the  ideal  4/3  power 
lav;,  i.e.. 


where  * is  the  coefficient  of  proportionality,  which  is  approximately 
2/3  -1 

0.01  enr'^  sec  . Similar  values  of  * for  average  conditions  existing 
in  the  sea  were  found  by  Inoue  (E.  .Inoue,  1952a),  Defant  (A.  Defant, 
1953), and  Hanzowa  (M.  Hanzcwa,  1953. 

Equation  (5),  obtained  experimentally  for  the  coefficient  of  tur- 
bulent exchange,  did  not  appear  unexpectedly,  because  as  early  as  1941 
A.  M.  Obukov  acquired,  for  the  coefficient  of  turbulent  exchange  from 
dimensional  considerations,  the  "4/3  power"  law,  which  is  based  cn  a 
theory  for  the  onset  of  turbulent  flow  proposed  by  A.  N.  Kolmogorov. 
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In  1960  Ichye  and  Olson  (T.  Ichye,  F.  C.  W.  Olson,  I960)  made  an 
attempt  to  extend  the  "4/3  power"  law  to  processes  with  a scale  of  up 
to  1000  km.  At  the  same  time,  Joseph  and  Sendner  (J.  Joseph,  H.  Sendner,\ 
1958)  and  similarly  Schonfeld  (J.  Schonfeld,  1961)  hypothesized  that 
for  such  macroscopic  processes  in  the  ocean  the  "4/3  power"  law  loses 
strength,  and  in  that  case,  the  coefficient  of  exchange  is  best  related 
to  l , i.e.,  a linear  dependence  exists.  That  question,  obviously, 
still  requires  further  investigation;  however,  even  now,  with  the 
dependence  of  A(l)  on  the  scale  of  phenomenon  established,  there 
is  the  possibility  of  being  more  correct  than  formerly  in  the  selection 
of  values  for  the  exchange  coefficient  for  the  solution  of  various 
problems  about  oceanic  currents  and  processes  of  exchange . 

Other  recently  developed  investigations  of  turbulence  in  the  sea 
are  based  on  the  statistical  and  spectral  theory  of  turbulence,  and 
their  principal  attention  is  devoted  not  to  investigations  of  average 
fields  of  various  substances  and  speeds  in  the  ocean,  but  to  the  study 
of  instantaneous,  fluctuating  values  of  these  quantities.  This  approach 
is  based  upon  the  representation  of  turbulent  flows  as  the  superposition 
of  many  scales  of  fluctuations,  i.e.,  the  generators  of  the  distinct 
turbulent  spectrum.  In  the  semi -empirical  theories,  turbulence  is  not 
actually  investigated;  thus  the  nature  of  turbulent  exchange  coefficients 
remains  vague  to  a significant  degree. 

The  distinct  statistical  characteristics  of  turbulence  have  been 
acquired  to  a sufficiently  broad  extent  in  recent  years  by  investiga- 
tions of  turbulent  flews  in  aerodynamic  tunnels  and  other  laboratory 
equipment,  and  also  in  layers  of  the  earth' 3 atmosphere.  Unfortunately, 
in  investigations  of  oceanic  turbulence,  the  statistical  methods  of 
investigation  are  quite  poorly  developed.  But,  before  dwelling  on  the 
works  bearing  on  that  question,  the  requirements  for  instrumentation 
for  measuring  fluctuations  of  quantities  in  a field  of  oceanic  turbu- 
lence will  be  briefly  described. 

If  the  average  speeds  of  currents  have  small  values  and  if  the 
periods  of  fluctuations  are  likewise  not  greaa,  for  measuring  the 
fluctuations  of  speeds,  it  is  obviously  necessary  to  have  a quite 
precise  and  lew  -inertia  apparatus.  This  is  particularly  germane  to 
the  study  of  vertical  turbulence,  which  is  a case  of  extremely  small 
average  speed,  and  wherein  the  scale  of  the  phenomenon  is  also  com- 
paratively small. 
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The  case  is  somewhat  different  for  investigations  of  horizontal 
turbulence  in  the  sea.  The  huge  horizontal  dimensions  of  the  seas  and 
oceans  provide  extremely  large-scale  horizontal  processes.  As  for  the 
average  horizontal  speeds,  so  the  fluctuations  of  these  speeds  have 
incomparably  larger  w*  nitudes  than  the  analogous  vertical  characteristics. 
For  the  study  of  horizontal  turbulence  in  the  sea,  the  requirements 
for  precision  and  "inertialness"  of  the  apparatus  employed  ought  to 
differ  from  those  used  for  measuring  the  characteristics  of  vertical 
turbulence.  If  in  the  latter  case  a sufficient  representation  of 
average  values  and  fluctuations  of  speeds  is  acquired  by  records  from 
low -inertia  devices  over  a period  of  a few  minutes, then  for  studies  of 
horizontal  turbulence,  in  order  to  be  able  to  more  or  less  envelop 
large-scale  changes  in  horizontal  speed,  it  is  necessazy  to  have  lengthy 
series  of  observations  (on  the  order  of  several  days) . But  the  inertial 
and  precision  requirements  of  these  instruments  may  be  significantly 
relaxed.  So  devices  which  record  speeds  averaged  for  2-3  minutes  as 
"instantaneous  speeds"  are  recognized  as  fully  satisfactory  for  the 
investigation  of  horizontal  turbulence  in  the  sea,  if,  of  course,  one 
is  not  interested  in  the  especially  subtle  microstructure  of  fluctua- 
tions, which  give  only  insignificant  contributions  to  the  general 
energy  of  horizontal  turbulence  in  the  sea.  For  studies  of  the  vertical 
turbulence,  such  devices  are  useless  and  a low-inertia  device  of  one 
type  or  another  must  be  designed. 

A well-known  success  in  this  direction  v/as  attained  in  the  Depart- 
ment of  Marine  Physics  of  the  Moscow  Government  University  under  the 
direction  of  A.  G.  Kolesnikov  (1958,  1959c).  With  the  aid  of  the  low- 
inertia  devices  (turbulence  meters),  measurements  of  fluctuations  of 
speeds  and  temperatures  were  conducted  in  the  cold-water  layer  in  Lake 
Baikal  and  at  the  stations  of  "North  Pole-4" . This  work  was  also 
conducted  in  the  Indian  and  Atlantic  Oceans  and  the  Black  Sea.  As  a 
result  of  the  analysis  achieved,  the  papers  of  A.  G.  Kolesnikov  and 
his  co-workers (1959a,  b,  I960;  N.  A.  Panteleev,  1959)  acquired  curves 
for  the  distribution  according  to  depth  (s)  of  turbulent  intensity, 
coefficients  of  turbulent  viscosity  >'(z),  and  thermal  conductivity 
k (z) . It  was  shown  that,  in  conditions  of  stratification,  the  average 
magnitude  of  v(z)  is  many  times  larger  than  the  value  of  k (z),  so 
that  on  the  average  v(z)  = 50 -k  (z) . They  also  discovered  the 
dependence  of  the  coefficients  v(z)  and  k (z)  upon  the  stability 
criterion  for  layers  of  water;  however,  a precise,  analytical  form 
for  this  dependence  was  not  achieved.  A.  G.  Kolesnikov  also  succeeded 
in  confirming,  in  sea  conditions,  the  well-known  law  of  A.  N.  Kolmogorov 
concerning  the  change  in  speeds  of  turbulent  fluid  motion  with  distance 
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according  to  a power  law,  with  a demonstrated  power  of  l/3.  It  is  to 
be  noted  that  all  these  results  v/ere  acquired  in  the  surface  layers 
of  the  sea,  while  the  continuous  determination  of  corresponding 
quantities  did  not  exceed  several  minutes. 


A device  for  measuring  speed  fluctuations,  based  on  the  electro- 
magnetic method,  was  applied  by  Bowden  and  Fairbaim  (K.  F.  Bowden, 

L.  A.  Fairbaim,  1952)  to  the  investigation  of  turbulent  fluctuations 
of  the  horizontal  speed  w in  strong  tidal  currents  at  two  meters  above 
the  bottom.  Calculations  of  autocorrelation  and  spectral  functions  for 
the  u component  of  speed  were  produced,  and  estimates  of  the  scale  of 
turbulence  were  also  made.  In  other  work,  Bowden  and  Fairbaim  (K.  F. 
Bowden,  L.  A.  Fairbaim,  1956)  described  experiments  of  simultaneous 
measurements  of  speed  fluctuations  « and  to  , that  made  it  possible 

for  them  to  compute  the  Reynolds'  stress:  p/Tv  . However,  the 

length  of  the  records  of  speeds  ir.  the  experiments  of  Bowden  and 
Fairbaim  did  not  exceed  5-10  minutes  because  they  acquired  values 
related  to  the  microstructure  rather  than  the  large-scale  structure; 
but-  this  flow,  as  already  shown,  has  a highly  specialized  form. 


Turbulently  pulsating  speeds  in  current  flows  near  Cape  Cannon 
were  measured  by  Nan'niti  (T.  Nan'niti,  1956)  with  the  help  of  his 
invention  of  a photoelectric  instrument  for  measuring  currents.  It 
was  revealed  that  the  ratio  of  average  amplitudes  of  pulsation  to  -the 
magnitude  of  average  speed  (the  intensity  of  turbulence)  was  . maximum 
near  the  bottom  and  minimum  in  the  middle  of  the  flow.  Computation 
of  the  autocorrelation  function  R(t)  showed  that  it  maintains  an 
absolute  law  for  various  series  of  observations  and  that  periods 
averaged  from  10  - 20  seconds  to  5 - 10  Canutes.  The  theoretical 

expression  for  the  autocorrelation:  1 - R (t)  ^ t2/3,  follows  from 
similarity  theory,  which  is  rarely  fulfilled  by  data  of  the  experiments 
described.  For  this  reason,  in  the  opinion  of  Nan'niti,  it  appeared 
that  the  duration  of  the  observations  and  period  of  averaging  were 
insufficient,  and,  likewise,  that  the  sensitivity  of  tne  device  was  low 


A low-inertia  device  of  the  thermal  anemometer  type  for  the  measure 
ment  of  turbulently  pulsating  speeds  was  applied  by  Grant,  Stewart, 
and  Moillet  (H.  L.  Grant,  R.  V/.  Stewart,  A.  Moillet,  I960)  in  a strong 
tidal  current  in  Discovery  Straits  on  the  coast  of  British  Columbia. 

The  device  measured  turbulent  fluctuations  with  frequencies  up  to  one 
kilocycle,  so  that  for  average  speeds  of  flew  of  100  err/ sec,  there 
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corresponded  a minimum  dimension  of  1 millimeter  for  the  decay  of 
turbulent  motion.  The  spectral  function  for  the  investigated  turbulent 
flows  was  computed,  and  it  was  shown  that  these  spectra  correspond  to 
the  theoretical  expressions  of  A.  N.  Kolmogorov  and  Heisenberg. 

In  addition  to  the  enumerated  investigations  of  the  microstructure 
of  the  speed  field  in  oceanic  turbulence,  there  is  a series  of  works 
dedicated  to  the  study  of  pulsating  water  temperatures  (L.  Liberman, 

1951;  E.  Inoue,  1952;  W.  W.  Englich,  1953;  T.  han'niti,  1957;  N.  V. 
IKontobqyzeva,  1958;  and  others).  These  works  revealed  a series  of 
laws  of  temperature  pulsations  in  the  sea  (especially  in  the  fozm  of 
the  autocorrelation  function)*,  these  laws  are  important  for  the  specific 
investigation  of  turbulent  heat  conduction  and  the  formation  of  tem- 
perature fields  in  the  sea. 

A significantly  smaller  effort  has  been  dedicated  to  the  investiga- 
tion of  horizontal  macroturbulence  in  the  sea  through  the  aid  of 
observations  over  horizontal  speeds  of  currents  by  ordinary  oceanographic 
rotor  devices.  The  first  works  in  this  province,  as  also,  in  general, 
in  all  regions  of  the  application  of  statistical  methods  to  the  study 
of  oceanic  turbulence,  appeared  in  the  investigations  of  V.  B.  Shtokman, 
accomplished  in  the  years  1940  - 1941.  In  the  first  v/ork  (1940),  V.  B. 
Shtokman  determined  coefficients  of  horizontal  turbulent  friction  in 
the  Caspian  Sea  using  data  of  current  speeds  measured  by  a hydrological 
rotor;  10 -minute  averages  were  formed  every  two  hours  in  the  course  of 
23  hours.  Coefficients  of  horizontal  turbulent  friction  were  computed 
by  the  formulae  of  Ertel,  which  were  generalized  expressions  of  the 
semi -empirical  theory  of  Prandtl  in  tensor  form.  By  a long  path  of 
synthesis,  V.  B.  Shtokman  determined  the  path  of  integration  for 
pulsating  speeds  in  the  limit  of  two  successive  changes  of  its  law. 

The  value  of  the  coefficient  of  horizontal  turbulent  exchange  appeared 

to  lie  within  the  bounds  of  1.2  to  4.9  X 106  cm2/sec,  with  an  abrupt 
change  of  value  for  anisotropic  exchange  relative  to  the  direction  of 
the  coast  line. 


In  his  work  of  1941,  V.  B.  Shtokman  used  the  following  formula 
for  the  determination  of  tne  horizontal  exchange  coefficient,  which 
follows  from  the  averaged  ecuations  of  fluid  motion: 
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With  this  formula,  V.  B.  Shtokman  found  A„  = 55  X 10^  cmVsec  and 

A„  =7  X 10^  cm2/sec  (the  x axis  was  directed  along  the  coast 
line) . 

The  simultaneous  measurement  of  current  speeds  with  several  rotors 
at  a series  of  points  in  the  vertical  enabled  V,  B.  Shotkman  to  compute 

the  coefficient  of  RH(*)  ( z ■ — the  vertical  coordinate!  As 

a result,  he  obtained  the  characteristic  curve  of  correlation,  which 
affirms  that  this  law  is  approximately  preserved  to  a depth  of  25  meters 
(the  sea  depth  was  35  meters  at  this  place) . 

Formula  (6)  is  the  same  as  that  used  by  Stomnel  (H.  Stommel,  1955) 
for  the  computation  of  the  coefficient  of  horizontal  turbulent  exchange 
in  the  Florida  Straits  by  data  of  speeds  measured  at  15 -minute  intervals 
from  an  anchored  vessel  in  1885.  Analogous  computations  were  produced 
for  the  region  of  the  Kuroshio  Current  (T.  Ichye,  1957),  using  a survey 
of  speeds  produced  by  the  aid  of  an  EMIT  (GEK) . 

Works,  repeating  to  a significant  degree  the  statements  of  the 
earlier  investigations  of  V.  B.  Shtokman,  were  acquired  for  the  Black 
Sea  by  A.  N.  Gezentzvey  (1961)  and  V.  V.  Hlopov.  Moreover,  A.  N. 
Gezentsvey  (1959)  examined  the  question  about  the  period  of  averaging 
fluctuating  speeds  on  the  acquired  values  of  exchange  coefficients. 

That  short  enumeration  of  work  is  limited,  unfortunately,  to  a 
list  of  investigations  devoted  to  the  study  of  the  characteristics  of 
horizontal  turbulence  in  the  sea  on  the  basis  of  the  measurement  of 
pulsating  speeds.  Determination  of  these  statistical  characteristics 
of  horizontal  macroturbulence  in  the  sea,  with  the  exception  of  the 
coefficients  of  exchange  and  some  magnitudes  computed  by  V.  B.  Shtokman, 
are  virtually  unrealized.  For  that  reason,  such  work  was  undertaken 
by  the  author  in  the  Dynamics  of  the  Sea  Laboratory,  Institute  of 
Oceanology,  AN  CCCP*  (R.  V.  Ozmidov,  1962). 

In  the  capacity  of  initial  material  for  the  determination  of 
statistical  characteristics  of  horizontal  macroturbulence,  we  used 
"data  which  was  measured  by  current  speed  rotors  for  many  days  at 
buoy  stations  and  t alien  in  the  Black  Sea  in  1956.  Buoy  stations  with 


^Translator's  Note:  Academy  of  Sciences,  USSR. 
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two  rotors  of  the  BFV-2  system  of  Alexseyev  at  20  and  100  meter  horizons 
were  deployed  at  a point  with  coordinates  42°  50'  N and  40°  25'  E and 
maintained  from  13  to  27  August . The  interval  between  the  sampling 
of  rotor  speeds  was  established  at  20  minutes.  Altogether,  the  tapes 
of  speed  records  consist  of  1008  imprints  of  vector  values  of  horizontal 
current  speeds  for  the  20-meter  horizon  and  1009  imprints  for  the  100- 
meter  horizon  of  observations. 


Ey  analysis  of  the  acquired  tapes  of  current  speeds,  it  was  not 
difficult  to  display  the  investigated  speed  field  in  the  presence  of 
irregular  oscillations  of  very  differing  magnitudes,  beginning  from 
Rtnai i scales,  which  provoked  the  analysis  of  successive  values  of 
speed,  and  finally  the  very  largest  scale  with  periods  commensurable 
with  the  length  of  the  entire  interval  of  observation.  These  large- 
scale  changes  of  speed  may,  of  course,  be  viewed  as  "regularities"  and 
may  be  studied  for  the  analysis  of  reasons  for  the  origin  of  the  flow 
and  its  variability  (for  example,  under  the  action  of  a changing  wind 
field) . However,  there  is  the  possibility  of  other  points  of  view,  in 
accordance  with  which  changes  of  vector  speeds  of  all  scales  may  be 
observed  as  "random1'  pulsations,  subject  only  to  investigation  ty 
statistical  methods.  The  first  such  view  of  large-scale  processes 
in  the  atmosphere  was  apparently  developed  by  Defont  (A.  Defant,  1921); 
he  suggested  viewing  even  such  large  foimiations  in  the  atmosphere  as 
cyclones  as  elements,  of  atmospheric  macro turbulence . A similar  point 
of  view  was  first,  and  also  quite  fruitfully,  developed  for  horizontal 
turbulence  in  the  sea  by  V.  B.  Shtokman  in  the  previously  cited  work. 
For  the  computation  of  statistical  characteristics  of  horizontal  tur- 
bulence in  the  Black  Sea,  we  proceed  to  those  representations  which 
view  all  speed  oscillations  for  all  dates  of  observation  as  random 
turbulent  pulsations  of  various  scales. 

In  that  case,  the  speed  component  of  currents  in  the  investigated 
flow  is  naturally  represented  in  accordance  with  formula  (l)  in  the 
f oitji  of  a sum  of  averaged  and  pulsating  speeds. 

Computation  of  the  average  speeds  fi  and  v for  all  periods  of 
observation  (14  days)  gave  the  following  magnitudes  (with  an  accuracy 

of  1 cn/sec) : u = 0,  0=8  cnv/sec  for  the  20-meter  horizon;  and 

« 0,  0=7  cn/sec  for  the  100 -meter  horizon  of  observation. 

Since  the  average  value  of  fluctuating  currents  and  V are  by 

their  definition  equal  to  zero,  as  a measure  of  these  amplitudes, 
there  were  computed  mean  square  values  of  the  fluctuating  deflections. 
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Computation  of  these  magnitudes  from  all  the  fluctuating  speed  values 
possessed  for  the  20 -meter  horizon  gave  the  following  figures: 

= 16.9  cn/sec  and  =12.8  cn/sec. 

Corresponding  magnitudes  for  the  100-meter  horizons  of  observation 

showed  the  following:  v(«')2  = 12.6  cm/ sec  and  - 

7.6  cn/sec.  Dividing  the  magnitudes  of  the  moan  square  value  of  speed 
fluctuations  by  the  average  speed  of  flow,  we  acquire  an  important 
statistical  characteristic  of  turbulence,  i.e.,  its  intensity  (longi- 
tudinal and  transverse):  ^ = 2.1,  I,  = 1.6  for  the  20 -meter 

horizon  and  /.=1.8,  /.  =1.1  for  the  100-meter  horizon. 

These  figures  show  that  the  intensity  of  turbulence  in  the  period  of 
observation  was  quite  high,  and  this  turbulence  was  not  completely 
isotropic,  since  its  longitudinal  intensity  v/as  less  than  its  transverse. 
As  anticipated  the  intensity  of  turbulence  decays  with  depth  — at  the 
100  -meter  horizon  it  is  significantly  belcw  the  intensity  on  the  upper 
or  20-meter  horizon. 

The  four  different  combinations  of  component  products  of  fluctuat- 
ing speeds  form  a symmetrical  tensor  of  correlation  moments,  the 
significant  components  of  which,  «'2  and  l/2  , were  adduced 

formerly;  then  the  components  n'tf=v'tS  were  computed  and 

gave  the  values  89  cm/sec  and  68  cn/sec*  for  the  20-  and  100-meter 
horizons,  respectively.  The  fact  that  the  non -diagonal  members  of 
the  tensor  of  the  correlation  moments  turned  out  to  be  different  from 
zero  again  confirms  the  anisotropy  of  the  turbulent  flow. 

As  a measure  of  the  statistical  tie  between  data  of  pulsating 
magnitudes  ( «'  and  v'  in  our  case),  the  so-called  correlation 
coefficient  R is  ordinarily  used;  it  is  formed  from  the  corresponding 
correlation  moment  divided  by  the  mean  square  deviation  of  the  observed 
pulsating  magnitudes.  Computation  of  correlation  coefficients  gave 
for  the  20-meter  horizon  R =-0.41  and  for  100-meter  horizon  R = -0.72. 

It  follows  that  the  pulsating  components  of  speed  were  more  random 
(chaotic)  on  the  upper  horizon  than  on  the  lower,  where  they  may  have 
.a  closer  relation  with  one  another. 


^Translator's  Note:  It  is  believed  that  the  units  intended  were 
(cn/sec) 2 vice  cn/sec. 
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In  addition  to  the  previously  mentioned  correlation  moments,  cor- 
relation moments  v/ere  computed  between  one  and  another  of  these  compon- 
ents of  pulsating  speeds,  measured  at  one  point,  but  at  a different 
moment  of  time: 


B.(0  -✓(<•>  (/.+<) 


(7) 


where  t„  is  some  moment  of  time  taken  as  the  initial  time,  and 
may  have  any  value. 

Since  the  initial  material  presented  was  not  a continuous  record 
of  speed  but  a series  of  discrete  values,  the  computed  correlation 
moments  (7)  were  not  produced  by  integration  in  time  but  by  means  of  a 
summation  over  a series  of  individual  samples  of  the  flow,  sampling 
from  seme  beginning  reading  through  a series  of  successive  moments  of 
time  on  the  tape  record  of  speeds.  The  number  of  the  many -valued 

functions  B.(t)  or  B,(t)  necessary  for  an  average  must  be 
foreseen  in  advance  to  avoid  difficulties,  and  for  us  it  was  selected 
as  equal  to  100,  i.e.,  each  value  of  correlation  moment  was  acquired 
as  a result  of  averaging  over  100  individual  values  of  the  product 

with  the  form  w't/,)  •«'(/,+/)  • For  the  investigation 

of  the  question  about  tne  significance  of  the  number  of  individual 

values  of  the  function  B.(t),  used  in  the  average,  we  produced  a 

special  computation  of  parts  of  the  function  B„(t ) for  different 

multiples  of  the  sampling  period.  That  computation  showed  that,  when 
the  toial  lag  period  ( *fc  ) equalled  from  50-75  times  the  averaging 
period,  the  function.  B.  CO  was  fully  developed  and  did  not  change 

by  a further  increase  of  the  lag  interval. 

For  the  clearest  graphs  of  the  computed  correlation  moment  func- 
tions, all  values  of  these  functions  were  normalized  to  unity, i.e., 
divided  by  the  product  of  the  mean  square  magnitudes  of  corresponding 
fluctuating  speeds  in  initial  and  in  sampled  moments  of  time: 
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The  functions  /?n(/)  and  R. (t)  were  computed  over  100 

samples  for  all  periods  of  observation  and  for  the  20-  and  100-meter 
horizons.  Graphs  thus  acquired  shot/  that  pulsating  speeds  are  correlated 
with  themselves  (autocorrelation)  greatly  over  lengthy  intervals  of  time. 
Indeed,  the  curves  for  /?.(/)  and  RAD  for  the  20 -meter  horizon 

approach  zero  as  a limit  for  the  first  time  at  t = 135  hours  and 
t = 81.5  hours,  respectively.  For  the  100-meter  horizon,  the  positive 
law  of  functions  RAD  and  R.(t)  is  preserved  even  further: 

141.5  and  150  hours.  But  with  the  expiration  of  such  a large  interval 
of  time  the  absolute  values  of  the  functions  RAD  and  RAD 

do  not  be cane  small  but  become  large  negative  values,  and  only  at  the 
end  of  the  observation  interval  when  the  argument  / 12  days  do  the 

functions  /?„(/)  and  RAD  again  approach  zero.  It  is  of  inter- 
est to  note  that,  for  the  20-meter  horizon,  the  first  point  after  the 
initial  point  of  functions  and  R RAD  shows  a fall  in 

correlation  of  26  percent  for  the  « component  of  speed  and  of  17 
percent  for  the  v component.  Likewise,  for  the  100 -meter  horizon 
the  drop  consists  of  7 and  2 percent,  respectively.  That  fact  states 
that,  in  the  upper  layer  of  the  sea,  a comparatively  large  role  is 
played  by  small  period  ( "t  < 20  minute)  fluctuations,  but  fluctuations 

for  <>  20  also  influence  the  autocorrelation  relationships  of  the 
given  speed  components.  At  the  greatest  of  these  depths,  such  3hort 
period  speed  fluctuations  are  developed  quite  weakly. 

In  addition,  for  the  basic  oscillations  with  a period  close  to 
the  periods  observed  on  all  graphs  of  functions  R.  (0 1 and  RAD 
for  the  20-meter  horizon,  it  is  not  difficult  to  display  (especially 
in  the  middle  part  of  the  graphs)  oscillations  with  periods  near  17 
hours,  i.e.,  near  to  the  period  of  inertial  oscillations  v/hich  may 
exist  in  the  given  region  of  the  sea.  For  a clearer  elucidation  of 
these  oscillations,  we  canputed  the  functions  /?.(/)  and  /?,(/) 
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for  the  20 -meter  horizon  by  choosing  the  starting  moments  of  time  very 
close  to  the  beginning  of  the  rise  of  the  phenomenon  studied.  For 
that  purpose,  as  initial  moments  of  time,  marks  an  the  tape  recording 
numbered  389  - 488  were  taken.  In  this  manner,  inertial  oscillations 
became  obvious  and  the  autocorrelation  functions  yield  a highly  accurate 
view  of  the  oscillating  process  with  decaying  amplitude  (Fig.  1).  The 
period  of  this  oscillation  equals  16  hours  and  40  minutes,  i.e.,  it 
corresponds  exactly  to  the  inertial  period  for  the  latitude  of  the 
observation  points.  Consequently,  the  construction  of  such  statistical 
characteristics  of  flows  as  autocorrelation  functions  makes  it  possible 
to  discover  the  presence  of  orderly  pavements  of  one  or  another  period 
in  turbulent  flows.  Some  theoretical  aspects  of  the  question  were 
examined  by  describing  other  statistical  characteristics  of  turbulence, 
namely, the  so-called  structure  functions  which  also  make  it  possible 
to  discover  the  presence  of  periodic  oscillations  with  inertial  period. 

As  already  shown,  the  turbulent  flow  studied  did  not  appear  isotropic, 
as  might  have  been  expected;  thus,  clearly,  the  large-scale  vortex 
formations  in  the  sea  ought , by  one  means  or  another,  to  depend  on  the 
geometrical  characteristics  of  the  flow  (coastline  features,  bottom 
relief) . In  the  capacity  of  a statistical  characteristic  of  such 
anisotropic  flows  one  ordinarily  uses,  as  first  introduced  by  A.  N. 
Kolmogorov,  tire  structure  functions,  determined  by  the  following 
expression: 


DJt=[r,  (/„-/)  -‘'.(b)  I 


(9) 


where  j , k equal  1,  2,  or  3,  ana  V)  and  vt  are  any  of  the  speed 
components . 

If  now  it  is  assumed  that  the  interval  of  time  t does  not  exceed 

scale  r , less  than  which  the  vortices  may  be  considered  isotropic, 
then  the  only  two  non -aero  structure  functions  are  shown  to  be: 

DJt) =|  f<(/.+0— ntf.TP ; 

D,  (<) = li/(/#+0  — (10) 
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In  scale  intervals  less  than  r and  greater  than  ~u  , where  the 
action  is  determined  by  viscous  forces,  for  structure  functions 
D,(f)  and  D,(t)  , it  may  be  found  from  dimensional  considerations 

that: 


DAD  : DAO^cSi'UW* 


(11) 


where  ct  and  c«  are  some  universal  dimensionless  constants  on  the  order 

of  unity,  < is  the  rate  of  dissipation  of  turbulent  energy  in  the  flow 
and  the  product  Ut  has  the  sense  of  "distance1'  between  the  "points" 
at  which  the  measurements  are  produced. 

Computations  of  the  structure  functions  D,{t)  and  D.(D 

were  produced  for  both  components  of  speed  and  for  the  20-meter  as  well 
as  the  100-meter  horizon  with  the  use  of  100  samples  in  averaging. 

The  results  showed  that  the  structure  functions  of  the  speed  fields 
in  the  sea  satisfactorily  agreed  with  formulae  (11).  Moreover,  there 
are  superimposed  on  the  general  dependence  of  the  "2/3  lav;"  some  large- 
scale  oscillations,  e.g.,  the  periodic  oscillations  with  an  inertial 
period  of  16  hours  and  4 0 minutes,  as  occurred  for  the  autocorrelation 
function.  We  will  examine  that  question  in  more  detail.  We  assume 
that  in  the  sea,  side  by  side  with  an  isotropic  field  of  chaotic  turbu- 
lent speeds  *>(t)  , there  exist  regular  components  of  speed,  changing 

according  to  a periodic  law:  Aca&*t  ; that  is,  the  total  speed  of 

flow  would  be: 


K(/)=y4CO*orf  + t>(/) 


(12) 
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Constructing  the  structure  function  for  such  a speed  field,  we 


(The  term  2[A  cos  mjju+t)  - A co»  «f0]  • M/,+0  - »(/<,)) 
as  a resultant  property  of  averaging) . 


(sin2a»(T+0  —sin  2*0 


d.{t)  is  the  structure  function  of  the  field  of  random 


where 

turbulent 


^Translator's  Note:  The  original  paper  gave  "t  where  it  should  be  "t? 
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If  the  averaging  period  T is  a multiple  of  the  period  of  oscillation 
of  the  regular  speed  component  (which  occurs  in  our  case),  then  we 
acquire: 


Z>r(0=2i4,sinl-£-+</.(0 


(15) 


The  structure  function  </,(/)  for  the  turbulent  speed  field 

t>(0  ought  to  obey  the  "2/3  law".  Consequently,  we  have  the 
conclusion: 


•t  . 


D,{t)  =2/lJsin2—  +c,e2'31/2,3<1;3 


(16)* 


The  graph  of  the  structure  function  computed  by  this  method  is 
presented  in  Fig.  2 in  which  experimental  points  of  the  naturally 
occurring  structure  function  for  the  v component  of  speed  at  the  100- 
meter  horizon  are  also  plotted.  The  graph  of  function  A(/)i  is 

constructed  by  formula  (16)  for  the  choice  of  the  following  values 
as  its  constants:  A = 2.7  cn/oec,  c+wu*1*  •»  1.72  X 10 

cn^/sec^ . The  magnitude  of  A is  determined  by  the  amplitude  of 
inertial  oscillations  existing  at  the  time  of  observations  in  the 
region  of  investigation  in  the  Black  Sea.  The  acquired  value  for  the 
magnitude  of  A is  in  good  agreement  with  the  value  of  the  amplitude 
of  inertial  oscillation  calculated  by  A.  D.  Yampolski  (1960)  with  the 
data  of  these  observations,  but  by  using  a method  different  from  that 
just  described.  The  values  of  the  two  constants  in  formula  (16) 

, make  it  possible  to  produce  an  estimate  of  an  important  characteristic 
of  turbulence:  the  rate  of  dissipation  of  turbulent  energy  « 


^Translator's  Note:  The  2 ft  ^ sin"2  term  appeared  in  the 

2 

original  paper  but  should  be  2 A ^ sin*  / 2. 
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Considering  that,  the  mean  speed  u for  the  investigated  segments  of 
time  was  equal  to  17.0  cm/sec,  ana  further  assuming  that  £ = 1,  we 

find  the  value  of  1.32  X 10"4  cm^/sec^  for  * . Although  direct 

measurements  of  « at  a depth  of  100  meters  were  not  produced  at  this 
time,  the  order  of  magnitude  of  c acquired  by  us  agrees  with  the 
available  indirectly  estimated  value  of  « in  the  sea.  (R.  V.  Ozmidov, 
1960) . 

The  good  agreement  between  the  theoretical  curve  and  experimental 
data  for  the  structure  function  and,  likewise,  the  realistic  acquired 
values  for  the  magnitudes  of  A and  « once  again  show  the  fruitfulness 
of  the  application  of  statistical  methods  of  the  theoiy  of  turbulence 
for  investigations  of  the  various  dynamic  processes  existing  in  the 
sea.  And,  undoubtedly,  these  methods  will  find  a quite  broad  applica- 
tion in  oceanographic  practice. 


20  METER  HORIZON 


TIME  IN  HOURS 


20  METER  HORIZON 


TIME  IN  HOURS 

Pig.  1 • Graphs  of  autocorrelations  functions  fy(t)  and  F^(t)  for  the  20  meter  horizon  and  for  the 
selection  of  Initial  times  as  marks  numbered  389  - 488 
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Pig.  2 - The  natrually  occurring  structure  function  CV(t)  for  the  100  meter  horizon 

1.  Theoretical  curve  of  the  *2/3  law.  * 

2.  Graph  of  theoretical  function  (16). 

3.  Experimental  values.  USL  Tech.  Memo.  No.  913-21-64 
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APPLICATION  OF  THE  THEORY  OF  RANDOM  FUNCTIONS 
TO  THE  STUDY  OF  OCEAN  CURRENTS* 

ty’ 

B.  N.  Belyaev  and  V.  S.  Boldirev 


In  the  past  fev/  years,  among  specialists  studying  ocean  currents, 
and  among  navigators  using  the  data  of  these  currents,  much  attention 
has  been  devoted  to  the  variability  of  ocean  currents  as  one  of  their 
essential  properties.  The  majority  of  investigators  have  shown 
L 3,  6,  9**_/  that  vector  current  speeds  appear  as  functions  of  time 
and  that  the  concept  of  constant  currents,  existing  up  to  now  in  the 
specialized  literature,  is  highly  conditional.  One  of  the  authors 
L*o J made  even  such  a conclusion  as:  "The  variability  of  speeds  and 
direction  of  currents  is  the  most  essential  of  their  characteristics 
in  any  part  of  the  sea  and  at  any  time  of  the  year.  It  may  be  ve in- 
definitely said  that  no  currents  are  absolutely  constant  for  any  in- 
terval of  time,  and,  in  seme  regions  of  the  ocean,  the  variability  of 
currents  has  so  much  significance  that  currents  of  any  direction  are 
equally  probable." 

To  a certain  degree,  the  variability  of  currents  is  undoubtedly 
conditioned  by  the  influence  of  tide -generating  forces,  i.e.,  there 
is  a tidal  nature  co  the  structure  of  current  vectors.  However,  if 
one  considers  currents  in  non-tidal  seas,  or  if  one  excludes  tidal 
components,  those  currents  without  such  components  similarly  preserve 
a characteristic  temporal  variability.  Physically,  this  is  easily 
explained.  Actually,  currents  are  agitated  by  many  factors  acting  in 
common.  All  these  factors  are  quite  variable  v/ith  time,  such  that  it 
is  said  to  be  impossible  to  use  their  essential  features  as  the  prin- 
cipal means  to  determine  instantaneous  current  vectors.  Moreover,  one 
may  consider  that  the  speed  and  direction  of  currents  appear  as  random 
functions  of  time  in  general.  Qualitatively,  workers  hypothesize  that 


*The  article  appears  in  Oceanology . Vol.  3,  No.  6,  1963,  pp.  953-961. 
^Translator's  Note:  Original  gave  the  numeral  9 but  there  were  only 

6 sources  listed  in  bibliography. 
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it  is  possible  to  accept,  with  a sufficient  degree  of  precision, that 
currents  which  are  examined  for  such  intervals  of  time  that  seasonal 
changes  of  general  hydrometeorological  conditions  are  insignificant 
appear  as  stationary  random  functions  of  time, for  the  physical  condi- 
tions of  stationarity  are  fulfilled  in  that  case. 

On  the  basis  of  observations,  it  becomes  clear  that  it  is  permis- 
sible to  apply  the  concepts  and  apparatus  of  the  theory  of  random 
functions  to  the  study  of  ocean  currents;  more  exactly,  its  most 
developed  section,  correlation  theory, can  be  used. 

As  is  well  known,  randan  functions  are  called  stationary  in  the 
broad  sense  of  the  word  (in  the  sense  of  A.  Y.  Hinchina),  if  their 
mathematical  mean  and  variance  are  constants,  and  if  their  correlation 
function  depends  only  on  the  difference  of  their  arguments.  For  real 
randan  functions,  the  conditions  of  Hinchina  are  nob  exactly  realised 
as  a rule,  and  sane  approximations  are  made. 

As  is  well  known,  the  basic  characteristics  of  stationary  random 
functions  of  time  , l({)  appear  in  its  mathematical  expectation 

Ml{ty  and  its  correlation  function  B(x)  ( t is  the  lag 

time).  These  characteristics  are  defined  by  the  formulae: 


m = Af*(/)=  lim  UlWdt,  (1) 

r—oo  T j 


B(x)  = Ml(t)l(t  + x) 


lim 

T—ao 


T 


(2) 


or,  approximately, 


m =s  I £ &<»>(* A). 


( i • ) 


2b 
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5 M : « 77  S &U)  (feA  + T)  Vl)  (M).  (2 ' ) 

k=*  l 


Here  A = ^ > g<i>  is  a sample  of  the  process;  T is  the  sample  length. 

Knowledge  of  these  characteristics  allows  one  to  decide  a series 
of  investigative  and  practical  tasks.  In  the  present  article,  two  of 
these  will  he  examined:  (1)  prognosis  of  currents;  (2)  apportionment 
of  the  tidal  composition  of  total  currents. 

For  these  tasks  we  used  material  from  multi-day  buoy  stations 
taken  in  western  parts  of  the  Atlantic  Ocean  in  the  process  of  expedi- 
tionary work  in  1960.  Computations  of  six -day  graphic  projections  of 
vector  current  speeds  at  a meridian  and  parallel  were  made  in  order 
that  all  reduced  data  be  contained  in  the  vector  of  current  speeds, 
which  would  not  be  obvious  in  further  specialization.  Computations 
were  produced  on  an  electronic  correlator  with  magnetic  recording 
constructed  by  A.  N.  Shakova,  and  as  a result  20  correlation  functions 
were  obtained,  the  normalized  average  values  of  which  are  shown  in  the 
diagranft  These  functions  correspond  to  currents  measured  over  five 
days  at  separations  of  100  - 300  miles  and  at  levels  of  25,  50, and  100 
meters.  All  curves  acquired  were  good  approximations  of  the  expression: 


B (t)  = CeraW 


O) 


Here  C is  the  variance  ; a and  p are  constants . (Actual  values  of 
« and  p lie  in  intervals  from  0.15  to  0.60  per  hour  and  from  0.33 

to  0.70  per  hour,  respectively.  The  data  obtained  were  sufficient  to 
confirm  that  the  correlation  function  of  the  currents  in  similar  regions, 
at  similar  times  of  the  year,  and  at  similar  depths  always  has  the 
form  (».) 


^Diagram  appears  at  the  end  of  the  text. 
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PROGNOSTICATION  OF  CURRENTS 

From  what  has  been  stated  it  is  clear  that  none  of  the  texts, 
like  atlases,  charts,  etc.,  are  able  to  satisfactorily  answer  a question 
which  has  important  significance  for  navigators:  Y/hat  are  currents  at 
a certain  point  at  a particular  moment?  Texts  give  only  the  average 
value  of  currents — their  mathematical  mean  from  which  deviations  of 
actual  values  may  be  quite  great.  It  is  perfectly  clear  that  the 
question  may  be  fully  answered  only  by  the  instantaneous  measurement 
of  currents  and,  in  the  spaces  between  measurements,  there  arises  the 
task  of  prognosticating  currents. 

Prognosis  of  currents  may  be  realized  by  the  basic  method  of  linear 
extrapolation  of  random  functions.  This  method,  worked  out  by  A.  N. 
Kolmogorov,  N.  Vinerom,  and  others,  is  highly  reliable  and  successful 
from  the  point  of  view  of  utilization  of  the  appendix  to  the  work  of 
Yaglana  £4_/  , from  which  we  took  all  the  useful  mathematical  apparatus 
used  belav. 

As  is  well  known  /_  4_/,  the  correlation  function  (3)  is  related 
to  the  spectral  density: 


(4) 

/ (X)  = - , 

X«  -f  2 aX*  -f  b* 

where 

<7  = 0*  — P»;  b = Ya*  + P*;  D = 2—  (o*  + p»). 

ft 


By  (4)  we  find  the  so-called  spectral  characteristic  extrapola- 
tion ‘ the  following  conditions  ought  to  be  satisfied: 


a.  4», (1)  analytic  in  lower  half- plane  and  grows  more 

slowly  in  this  half --plane  than  some  power  of  jX|  as  | X I -»  * 
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t>.  i't>t(l)  = [cat-0>x(l))/(X) 


analytic  in  the  upper  half -plane  and,  as  IX  |— ► oo  , 

decays  in  this  half -plane  more  quickly  than  |X|— »— «t  e>0; 


c-  J l®T(l)|V(X)dl<oo. 


that: 


It  is  not  difficult  to  show  that  for  these  conditions  it  follows 


Ot  (X)  — Ak  + B, 


(5) 


where  constants  A and  B are  found  from  the  system: 


e‘Xt  - (X)  = 0, 

for  Xt  a 3 - r i'a  and  X,  = — (P  — ia). 


(6) 


After  calculations  and  transformations,  we  obtain: 


.X 


Ot  (X)  = i — r-°T  sin  (Jt  + e~a*  ^cos  pt  4-  sin  Pt  j 


(7) 


Values  for  the  randan  functions  related  to  <ht(X)  are 
extrapolated  by  the  expression: 


!«*»(/+ T).  J #«MI>x(X)dz(X). 


(8) 
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Substituting  (7)  into  (8)  and  having,  as  is  known  in  the  theoiy 
of  random  functions,  the  equalities: 


and 


(9) 


$ (0, 


(10) 


where  |(/)  is  the  value  of  the  random  function  at  the  instant  t , 
and  6' (0  is  its  derivative  at  that  same  moment,  we  finally  find: 


I<«  (t  + 1)  - e~at  t'“>  (0  + (cos  pr  + *p  sin  pr  j £<‘>  (/)} , (11) 


or 


V«v  + »)-«“•«  ^ V»  (t)  + 3 (t)  l (t).  (12) 
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If  M\(t)  * 0,  then,  before  starting  to  prognosticate, 

Mfc(f)  is  normalized  to  have  a zero  mean. 

Ey  formula  (12)  we  produced  the  prognosis  of  currents  at  3 points 
in  the  Atlantic  Ocean  at  1,  2,  3,  and  4 hours  ahead;  for  each  sample 
there  were  calculated  five  prognosticated  values.  For  all  there  were 
acquired  240  extrapolated  values,  each  of  which  corresponded  to  an 
actual  value  known  from  measurements. 

For  example:  the  predicted  currents  for  one  of  these  stations 
are  shown  in  Table  1. 

The  following  differences  were  formed  judging  the  effectiveness 
of  prognosis: 


A1  = S(/+t)-S(/  + t), 
A ,-*(/)-£(<  + *)» 


W+t)  is  the  prognosticated  value  of  the  current; 

is  the  actual  value  of  the  current;  £(/)  is  the 
actual  current  value  at  the  initial  time. 

Later  and  A2  were  expressed  in  terms  of  or  , the  mean 

square  deviation  of  the  data,  when  we  assume  the  instantaneous  value 
of  current  to  be  equal  to  its  average  value,  and  by  the  formula: 


where 


■(tr- »] 
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There  were  computed  the  error  in  prognosis  and  the  error  for  the 
case  of  the  current  at  moment  i<+*  having  the  value  it  had  had  at 
the  beginning  moment  / . 

Results  of  these  computations  are  shown  in  Table  2. 

Moreover,  if  at  the  moment  <+*  there  is  at  our  disposal  the 
possibility  of  selection  from  three  values  of  current:  the  averaged 
Mi(t)  , the  measured  £(/)  , and  the  prognosticated 

|(f+t)  , the  most  accurate  of  these  is  considered  to  be  the 

prognosticated  value  (remembering  that  we  always  understand  the  value 
of  the  vector  projection  of  current  speed  on  a meridian  or  parallel 
by  the  words  "current  value".) 

Of  course,  so  far,  this  method  is  to  be  applied  to  currents  whose 
tidal  constituents  are  sufficiently  small  or  have  been  excluded  from 
the  total  current.  Our  further  dissertation  will  be  dedicated  to  find- 
ing ways  which  may  be  used  to  exclude  the  tidal  constituents. 

ISOLATION  OF  TIDAL  COMPONENTS  IN  TOTAL  CURRENTS 

In  general,  we  represent  currents  in  the  form  of  the  sum  of  three 
functions : 


£(0  = Af&(')  + tj(0  + «P(0. 


(17) 


USL  Tech  Mono 
No.  913-27-64 


Not  diminishing  the  generality  of  our  development,  we  may  allow: 


then 


6(0  “»l(0  + <p  (<): 


(18) 


the  function  <p(t)  may  be  placed  in  a Fourier  series  expansion: 


where 


9(0  = 2 Ak  sin  ( kat  -f-  8*). 


(19) 


(i)  — ■“  *i  Ai,  = Y a\  b\ ; tg  C*  = — ■ ; 

T bk 


a* 


and  bk  are  the  Euler  coefficients  of  the  Fourier  expan- 


sion of  the  function  <p  (/).  . In  this  expansion,  U would  represent 


Mm 


It  is  easy  to  show  that  as  a result  of  the  treatment  of  (17)  on 
a correlator,  or,  similarly,  inserting  (17)  into  formula  (2)  and  (2') 
we  acquire: 


°°  A* 


J3*.  (t)  =*  (t)-f  V --  cos  I; u>r 


* i 


(20) 


The  correlation  function  /3n(r'  of  the  randan  function 

ti(/)  tends  to  zero  as  t->  <*  (in  practice,  when  t = 15  to  20 

are 


hours,  ’ members  of  we  — cos  feat 


r ~ * 

",  *1 
v'  —l 


undamped,  and,  consequently, 


fli(t)  is  likewise  not  damped  as 


T-*  oo. 
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Moreover,  the  presence  of  an  undamped  periodic  component  in  the 
graph,  achieved  as  a result  of  treatment  on  a correlator  of  sample 
current,  appears  as  a necessary  and  sufficient  indication  of  the 
existence  of  a tidal  current  at  the  investigated  point. 

On  the  contrary,  the  absence  of  a periodic  undamped  component  in 
the  composition  of  Br  (ti  , . appears  as  certification  that  the  tidal 

current  is  absent  at  the  data  point  and  that  the  variability  in  the 
current  is  formed  only  by  random  causes. 

In  the  simplest  case,  when  a hodograph  of  vector  speed  of  tidal 
current  is  presented  as  an  ellipse  describing  the  vector  current  with 
constant  angular  speed  « (which  may  apparently  be  considered  real  in 
the  open  ocean  at  great  distances  from  the  coast  and  at  great  depths), 
it  is  not  difficult  to  shop;  that  the  undamped  component  in  equation 
(20)  has  the  form: 


COS  (OX.  (21) 


From  the  graph  of  (21),  it  is  easy  to  select  amplitude  a and 
frequency  u>  so  as  to  construct  the  function: 


(22) 

c sin  tot',  t'=t 

defining  the  tidal  current. 

Clearly,  the  amplitude  of  the  total  current,  |(/)  , is  always 
larger  than  the  amplitude  of  the  tidal  current  <r(0  . 

From  (21)  and  (22),  in  order  to  acquire  the  value  of  the  amplitudes  of 
the  tidal  current,  it  is  necessaiy  to  double  the  magnitude  of  amplitude 
; of  the  undamped  component  removed  from  the  graph  of  tS~  ’Tl  and  then 
extract  the  square  root.  The  result  is  the  exact  value  of  the  amplitude 
of  the  tidal  current  in  the  case  of  an  elliptic  hodograph  and  an  approxi- 
mation which  becomes  more  exact  the  closer  the  graph  of  current  is  to 
a sinusoidal  curve  in  all  other  cases. 
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In  treatments  of  our  realizations,  the  amp^’tude  of  the  projection 
of  the  speed  of  the  total  current  was  40  - 50  cn/sec,  when  the  amplitude 
of  the  tidal  component  was  not  in  excess  of  6 - 8 cnysec.  The  graph  of 

the  undamped  part  of  v B((x)  presented  a curve  which  was  quite  near 

to  a cosinusoid;  that  certified  the  nearness  of  the  graph  to  a sinusoidal 
curve,  which  in  its  turn  certified  the  nearness  of  the  hodograph's  form 
to  an  ellipse;  at  last,  the  acquired  value  of  the  amplitude  was  shown 
to  high  accuracy. 

If  we  examine  the  realization  as  a tidal  current  by  the  ordinary 
methods  of  harmonic  analysis  (for  example,  by  the  method  of  Darwin), 
the  presence  of  a tidal  current  with  an  amplitude  up  to  40  - 50  era/ sec 
would  be  determined.  Curiously,  the  same  magnitude  for  the  tidal 
current  would  be  substantiated  by  theoretical  computations  based  on  the 
determined  hydrodynamic  levels  in  the  ocean.  For  an  example,  we  go  to 
an  article  by  Boris  Z 1 _7.  In  that  article  there  is  a chart  of  the 
North  Atlantic  an  which  are  shown  the  computed  theoretical  amplitudes 
of  tidal  current  speeds.  For  the  region  in  which  the  observations 
employed  in  this  work  were  produced,  theoretical  computations  gave  the 
magnitude  of  the  amplitude  as  being  on  the  same  order  as  that  which  we 
used.  The  author  of  the  article  Z !_/  assumed  that  the  result  of 
computations  of  the  amplitude  does  not  coincide  with  calculations  from 
observations  because  of  a defect  or  insufficiency  in  the  theoretical 
conditions.  In  view  of  the  former  discussion,  it  may  be  considered 
that  the  demonstrated  discrepancy  appears  as  a consequence  not  of 
inadequate  theory,  but  of  inadequacies  in  the  method  of  processing  the 
experimental  data,  by  which  the  total  and  not  the  tidal  current  was 
subjected  to  harmonic  analysis.  V/e  noticed  that  errors  in  the  ampli- 
tudes of  the  projection  of  the  tidal  current  speed  inevitably  yielded 
errors  in  the  orientation  of  the  current  ellipse. 

Moreover,  harmonic  analysis  ought  to  be  preceded  by  the  processing 
of  the  realization  on  a correlator  for  the  purpose  of  estimating  the 
tidal  component.  From  the  preceding,  the  expansion  in  a Fourier  series 
of  the  undamped  component  ("tail")  of  the  correlation  function  (20)  at 
once  obtains  the  value  of  the  frequency  u and  the  amplitude  A*  of 
the  terms  in  the  Fourier  series.  It  is  more  difficult  to  find  the 
magnitude  of  (-)*,  , which  determines  the  position  of  the  function 

<rf ')  at  the  initial  moment.  That  step  may  be  determined  either 

by  the  method  of  successive  approximations,  or  by  the  fonnula  of  Wyrtki 

Za./. 


3C 
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The  method  of  successive  approximations  in  the  case  of  an  ellip- 
tical hodograph  consists  of  the  graph  of  function  (22)  combined  with 
the  outcome  of  realization  (18)  in  order  to  determine,  by  sight,  the 
most  probable  position  of  curve  (22)  on  the  time  axis.  After  the 
generation  of  the  difference  function: 


ih  (0-6(0 sin 


which  is  again  generated  an  the  correlator,  and,  if  the  combination 
were  successfully  accomplished,  there  should  result  the  acquired  cor- 
relation function  without  an  undamped  periodic  component.  In  cases  of 
adverse  combinations,  this  computation  is  repeated  for  some  displacement 
of  curve  (22)  along  the  time  axis;  after  that,  the  computation  of  a 
different  curve  r\2(t)  , is  produced  repeatedly  on  the  correlator 

until  the  correlation  function  is  completely  deprived  of  its  undamped 
member.  This  final  version  of  the  function  makes  it  possible  to  find 
the  tidal  component  from  the  simple  equality: 


<p  (0  = S(0 


which  is  then  subjected  to  harmonic  analysis. 

If  there  is  a Fourier  series  representation  (e.g.,  for  a very 
complex  hodograph),  the  method  of  successive  approximations  adapts 
itself  to  a full  analog,  the  only  difference  being  that  each  time  a 
first  harmonic  A\cnswl'  of  the  Fourier  series  is  deduced  from 

the  decomposed,  undamped  periodic  components  of  the  correlation  func- 
tions, which  are  acquired  after  n processings  an  the  correlator  of 
n- ! difference  functions . 

V.'yrtki 1 s formula  permits  the  immediate  determination  of  the  angle 
• , which  locates  the  position  of  the  corresponding  Fourier  series 

(or  curve  (22))  an  the  time  axis.  Therefore,  in  general,  the  number 
of  processings  of  the  curves  on  the  correlator  is  decreased,  but  the 
number  of  curve  tracings  is  necessarily  increased  for  the  following: 
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for  each  harmonic.  The  angle  Qk  is  acquired  by  the  formula 


where  the  corresponding  mathematical  means  for  the  proceeeed  curves  are 
determined  cn  the  correlator. 


T IN  HOURS 


» 


Ua  Tech.  Memo.  No.  013-27-64 
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SCME  RESULTS  OF  INVESTIGATIONS  OF  TIDAL  CURRENTS 
IN  THE  OPEN  PARTS  OF  THE  PACIFIC  OCEAN  * 

by 

R.  R.  Belevich 

Far  Eastern.  Scientific  Investigation  Hydrometeorology  Institute 


For  the  explanation  of  the  character  and  the  strength  of  tidal 
currents  in  the  open  parts  of  the  Pacific  Ocean,  the  l/S  U SHQKALSKII, 
an  her  second  cruise  (November  1960  - January  1961),  at  a line  of 
station  points  in  the  northern  central  regions,  produced  moderately 
broad  studies;  five  days  of  cyclical  current  measurements  were  analyzed 
by  the  scheme  proposed  by  A.  I.  Dubaninin-1- . Investigators  employed 
electromagnetic  current  meters  (EMIT) . In  the  course  of  a day,  the 
ship  cruised  by  a plan  of  two  circuits  every  three  hours.  Daily  cycles 
were  accomplished  for  nearly  four  days. 

Harmonic  constants,  found  by  the  conjugate  method,  gave  satisfac- 
toiy  agreement  between  calculated  and  experimental  data  in  70  percent 
of  the  cases. 

It  was  found  that  there  was  better  agreement  for  waves  of  semi- 
diurnal period  than  diurnal. 


*The  article  appears  in  Oceanology.  Vol  2,  No.  2,  1962,  pp.  260-261. 
-*-Dubaninin,  A.  I.,  1955,  "Certain  Characteristics  of  Tides  by  Short 
Duration  Observations,"  Tr.  GOV.  Oceancm-.  Inst..  Issue  30  (42). 
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Tidal  Current  Ellipse  in  Open  Oceans  by  Observations 
a - Diurnal  Wave;  b - Semi-Diurnal  ’Wave 

Magnitudes  of  current  speed  were  acquired  with  the  EMIT;  the  abso- 
lute correction  "K"  was  not  introduced. 

Results  shaved  that  tidal  currents  obtained  significant  speeds 
(up  to  60  - 70  cn/sec)  on  the  surface. 

The  character  of  tidal  currents  is  irregular, diurnal.  The  ellipse 
of  daily  currents  is  nearly  circular,  as  expressed  by  the  vector  diagram. 

Besides  evidence  from  this  five-day  period,  obtained  under  diffi- 
cult celestial  navigation  conditions,  data  were  taken  for  one-day 
periods  in  different  regions:  in  the  zone  of  the  Northern  and  Southern 
Tradewind  Currents  and  in  the  zone  of  the  Equatorial  Countercurrent. 
Results  of  these  investigations,  as  in  the  five  days  listed  formerly, 
confirmed  the  existence  of  different  types  of  tidal  currents  in  the 
indicated  regions  of  the  Pacific.  Observers  also  confirmed  the  exis- 
tence of  constant  currents  present  in  the  regions  cited,  and  only  the 
Equatorial  Countercurrent  shaved  few  deviations  to  the  north  and  in 
greatest  speed  as  compared  with  average  data  for  the  season.  Simul- 
taneous observations  of  hydrological  and  hydrodynamical  elements  shewed 
significant  oscillations. 


CONCLUSIONS 

1.  Tidal  currents  in  the  open  ocean  are  of  enough  significance 
to  cause  the  measurement  of  steady  currents  to  be  approached  carefully. 
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2.  With  EMIT,  it  is  possible  to  measure  tidal  currents  in  open 
and  deep  waters  for  the  purpose  of  using  their  harmonic  constants  to 
distinguish  residual  (constant)  currents. 

3.  In  regions  where  significant  tidal  currents  exist,  there 
occur  similar,  significant  oscillations  of  hydrologic  and  hydrodynamic 
elements;  it  is  necessary  to  take  cognizance  of  them  in  different  types 
of  computations. 


Submitted  to  the  Editor 
4 September  1961 
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ON  THE  PERIODS  OF  INTERNAL  WAVES  IN  DEEP  INLAND  SEAS* 

toy 

0.  M.  Belyakov,  U.  M.  Belyakov 
Black  Sea  Department,  Naval  Hydrophysical  Institute 
Acadeny  of  Sciences,  USSR 


In  the  past  decade  there  has  grown  a collection  of  works  about  one 
of  the  least  studied  oceanological  problems,  i.e.,  the  investigation  of 
the  character,  distribution,  and  origin  of  so-called  "internal  waves" . 

The  hydrosphere  presents  itself  as  a medium  consisting  of  a series 
of  water  layers  with  differing  density.  In  such  a layered  structure  of 
water,  when  the  interfaces  between  layers  are  under  the  action  of  var- 
ious excitation  forces  (tidal  forces,  atmospheric  disturbances,  etc.,) 
internal  waves  are  generated  with  amplitudes  attaining  decameters  and 
sometimes  hundreds  of  meters.  Naturally,  such  "appreciable"  vertical 
changes  in  the  water  mass  play  a large  role  in  the  dynamics  of  the 
world  ocean.  The  significance  of  internal  waves  as  factors  contribut- 
ing to  water  variations  and  the  transmission  of  various  hydrological 
characteristics  to  deep  layers  is  very  great,  not  only  for  the  hydro- 
logical  regime  of  basins  but  for  their  biological  domains  as  well  by 
aiding  the  distribution  of  various  things  having  importance  to  living 
organisms.  Knowledge  of  the  magnitude  of  intemeil  wave  parameters 
and  of  the  laws  of  their  variability  also  has  value  for  submarine 
navigation,  for  hydroacoustic  prognosis,  for  the  needs  of  the  fishing 
industry,  and  for  other  purposes. 

In  various  parts  of  the  world  ocean  there  exist  internal  waves 
vith  different  characteristic  parameters  depending  upon  their  dimen- 
sions, the  configuration  of  the  coasts,  and,  especially,  temperature  and 
salinity  fields,  which  create  the  density  stratifications  that  are 
typical  for  each  water  basin. 


* Published  in  Oceanology.  Volume  III,  No.  6,  1963  (UDK  551.465,  pp.  962-9). 


USL  Tech  Memo 
No.  913-27-64 


► 

Lvs 


The  present  work  contains  the  analysis  of  germane  data,  acquired 
at  two  multi-day  stations  in  one  of  the  deep,  inland  seas  at  distances 
of  6 and  28  miles  from  the  shore  for  the  purpose  of  determining  the 
periods  of  internal  waves.  These  stations  were  occupied  in  the  spring 
for  periods  of  4 - 5 days  and  were  synchronozed  between  two  ships  which 
were  deployed  along  a track  normal  to  the  shore.  There  was  a two-hour 
interval  between  observations. 

For  the  analysis  of  our  data  acquired  at  these  stations,  we  studied 
time-varying  vertical  oscillations  of  density.  We  took  the  average 
values  of  density  at  standard  horizons  for  the  period  of  observation. 
Moreover,  time-varying  oscillations  of  these  isopycnals  were  also  re- 
lated to  standard  horizons.  Similarly,  the  representation  of  the 
studied  quantities  makes  it  possible  to  compare  the  parameters  of 
internal  waves  at  one  or  the  other  of  these  horizons  in  various  regions 
of  the  sea  in  various  seasons. 

Deep  inland  basins  possess  special  conditions  for  the  generation 
of  internal  waves.  In  observed  internal  waves  in  the  sea,  tidal  forces 
do  not  appear  to  be  the  principal  factor  causing  oscillations  of  water 
particles.  The  greatest  amplitude  of  tides  here  is  an  the  order  of  10 
cm  (at  the  surface);  the  character  of  the  tides  is  irregular,  semidiurnal* 
Of  course,  the  existence  at  considerable  depths  of  significant  vertical 
oscillations  attaining  100  -150  meters,  as  seen  in  the  analysis  of  the 
data  of  the  expedition,  reflects  the  nan-tidal  origin  of  these 
fluctuations.  V.  V.  Schulaykin  and  A.  D.  Dubrovelcki  were  the  first 
to  show  the  seiche  origin  of  internal  waves  in  inland  seas. 

It  is  well  known  that  the  excitation  of  internal  v/aves  in  any  place  of 
the  sea  can  be  caused  by  small  impulses  due  to  changes  in  the  atmos- 
pheric pressure,  changes  in  the  wind  field,  seismic  shocks,  etc. 

Waves  caused  by  various  impulses  are  found  in  various  stages  of 
development  or  attenuation,  are  dispersed  in  all  possible  directions, 
and  are  superimposed  upon  the  unattenuated  periodic  oscillations  of 
tidal  character;  thus,  systems  of  internal  waves  are  generated  in  a 
given  basin.  It  was  quite  obvious  that,  even  for  a repeat  of  one  of 
these  synoptic  situations,  the  quantitative  state  of  the  spectra  of 
internal  waves  may  change  for  the  investigated  place.  At  the  same 
time,  it  seems  that  there  is  an  inherently  determined  combination  of 
periods  for  each  section  of  the  sea. 

Due  to  the  fact  that  internal  waves  by  their  nature  and  mechanisms 
develop  analogously  to  surface  waves,  we  recognized  the  possibility  of 
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using  one  of  the  methods  for  quantative  estimate  of  characteristics  of 
sea-surface  waves  for  the  determination  of  internal  wave  periods, 
namely,  the  method  of  the  theory  of  random  functions.  The  basic  char- 
acteristics of  random  functions  appear  as  the  average  value  (mathematical 
mean)  for  the  correlation  function.  It  is  to  be  noted  that  for  simple 
computations,  for  known  admissibility,  we  used  properties  of  what  are 
called  ” stationary”  randcm  functions.  In  this  case,  the  wave  function 
is  given  a zero  mean  value,  and,  for  the  characterization  of  the  random 
function,  it  is  sufficient  to  compute  its  correlation  (autocorrelation) 
function. 

The  correlation  function  R(x)  of  random  function  /(/)  ’ 

appears  as  a quantitative  measure  of  the  persistence  of  the  random 
function,  when  samples  have  been  separated  one  from  the  other  by  the 
magnitude  r ; it  is  approximately  determined  by  the  expression: 


T 

^ (t)  ~ + T)^ 


-r 


(1) 


The  value  of  the  correlation  function  for  x=  0 is  the  variance 
of  the  random  process. 

In  practice,  a small  number  of  values  for  the  analysis  of  surface 
as  well  as  internal  waves  are  available  to  perform  the  calculation  of 
their  correlation  function  and  are  related  to  the  Fourier  transform 
(in  agreement  with  the  well-known  theory  of  Hinchina)  for  the  so-called 
spectral  density  S(w)  . Spectral  density  and  correlation  functions 
are  related  to  one  another  by  their  Fourier  transforms  as  follows: 
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(2) 

O) 


4B 


USL  Tech  Memo 
No.  913-27-64 


where  w is  the  angular  frequency.  The  function  S(io). , or  spectral 
density,  appears  sis  the  frequency  characteristic  of  the  randan  function 
f(t),  determining  for  the  process  its  distribution  of  energy 
according  to  spectral  frequency. 

In  our  case  the  randan  function  f(t)  is  mixed  with  the  undamped 
periodic  fi(0  . The  correlation  function  of  the  sum  f{t)  + 

/,(/)  will  also  be  the  sum  , in  which  /?(t) 

would  deccy  to  zero  with  an  increase  in  r , and_  (t)_  will  be  periodic 
with  the  frequency  of  the  germane  function  Lzj  and  /_  4 In  the  case 
of  mixed  functions  formula  (2)  is  also  applied. 

A mechanical  correlator  was  used  for  the  acquisition  of  the  cor- 
relation functions,  since  numerical  determination  is  tied  to  very 
unwieldy  computations . By  polar  graphical  representation  of  our  time 
functions  by  two  indices  with  changing  r in  friction  mechanism  of  the 
correlator  we  created  the  continuous  summary  product  /(/) 
and  in  our  counters,  after  the  circular  curve,  we  acquired  a magnitude 
proportioned  to  the  integral  ^ f(t)-f(t  + r)dx,  where  T is  the 

duration  of  the  function  in  hours.  » The  approximate  value  of  the  cor- 
• relation  function  is  acquired,  finally,  if  we  multiply  the  value  of  the 

l 

coefficient  on  the  device  by  the  value  r—x  . The  integration 
- steps  were  taken  as  1 hour  40  minutes. 

Computation  of  ' S(ia)  was  produced  for  the  ordinate  of  density 
temporal  oscillations  relative  to  the  average  horizon  at  25,  50,  75, 

100,  and  150  meters  for  the  offshore  station  and  at  25,  100,  150,  200, 
300,  and  422  meters  for  the  inshore  station.  For  this  computation, 
the  following  approximate  formula  was  used: 

r 

5 (w) « — ^ R (t)  cos  axdx.  (4 ) 
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The  integration  was  produced  numerically  by  the  formula 


2 R„  cos  an  A/  + 


cosrnmA/ 


where  Rt  is  the  value  of  the  correlation  function  for  t=  0; 
is  the  value  of  the  correlation  function  over  equal  intervals  of 
and  m is  the  quantity  of  pieces  A t 


Since  « = — - , for  the  convenience  of  computations  of  S(a>) 

we  changed  the  value  of  T . Up  to  T = 22  hours,  AT  was  taken  to  be 
5 minutes;  for  T>  22  hours,  AT  v/as  taken  as  1 hour  40  minutes.  For 
the  precise  location  of  the  maximum  spectral  density,  the  intervals 
AT  v/ere  increased  two  and  even  four  times. 


From  the  relations  (2)  and  (4)  it  is  quite  obvious  that  the  exact 
computation  of  the  spectral  density  depends  basically  on  the  exact 
determination  of  the  correlation  function.  The  latter  clearly  depends 

on  the  maximum  value  of  t , the  length  of  integration  T athe  magnitude 
of  the  integration  step,  and  on  the  instrumentation  errors  of  the  apparatus. 
As  far  as  v/ork  on  the  mechanical  correlator  is  concerned,  the  two  operations 
provide  an  opportunity  for  very  large  errors  to  arise  from  calculations 
because  of  the  subjective  character  of  compression. 

Because  of  the  approximately  determined  correlation  functions, 
errors  in  observation,  etc.,!  S(<o)may  have  a negative  value, 

2 n 

which  it  should  not,  for  some  value  of  «=y.  . The  spectrogram 

in  Fig.  1 may  serve  as  an  example.  The  value  of  (<o)  is  negative 
at  T,  equals  5.8,  7.5,  16.6,  and  42.0  hours.  Clearly, for  the  absence  of  the 
enumerated  errors  at  some  of  these  frequencies,  the  ordinate  s( to) 
ought  to  be  either  equal  to  zero  or  nearly  equal  to  zero.  However, 
in  spite  of  these  deficiencies,  which  occur  in  c amputating  the  magnitude 
of  the  spectral  density  with  the  aid  of  the  method  summarized,  it  is 
possible  to  acquire,  from  the  spectra,  the  periods  of  the  complex 
oscillations  which  compose  the  waves  and  which  we  observe  in  nature. 
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The  accompanying  tables  present  all  the  significant  periods  we  1 
acquired  in  the  interval  of  6 hours  to  3 days  at  the  offshore  and 
inshore  stations.  (Tables  1 and  2). 

Internal  wave  periods  are  quite  well  grouped  according  to  the 
vertical.  Ten  isolated  groups  of  periods  are  distinguished  by 
characteristic  values.  The  short  periods  - 5.8,  6.7,  7.5,  and  8.5 
hours  - belong  to  Group  I.  They  were  present  at  both  stations  from 
25  to  422  meters  deep.  These  periods  compare  completely  with  the 
well-known  periods  for  seiches  of  the  sea  surface  in  this  region. 

Group  II  contains  the  components  with  periods  of  9 - 10  hours,  which 
are  also  of  seiche  origin. 

Group  III  is  of  great  interest  because  its  periods  conform  to 
periods  of  semidiurnal  tidal  character,  as  observed  in  the  interior 
of  the  seas.  In  spite  of  the  insignificant  ( compared  to  the  former) 
amplitudes  of  the  surface  tides,  they  were  later  observed  to  be 
influential  in  all  layers  of  water.  From  the  examination  of  the 
tables,  it  is  seen  that  waves  with  a semidiurnal  period  exist  at  all 
the  horizons  we  investigated. 

The  existence  of  waves  with  periods  of  Group  V is  similarly 
predicted.  They  are  oscillations  with  the  inertial  period.  In 
agreement  with  formula  ( T = ■ -1-2  *,  the  period  of  internal  inertial 

\ sin<p  / 

waves  for  the  average  latitude  equals  17  to  18  hours. 

Further,  from  the  analysis  of  the  tables  it  is  to  be  noted  that 
the  nature  of  the  inshore  and  offshore  stations  characterizes  their 
groups  of  waves.  Thus,  at  tne  offshore  station  there  exist  components 
with  24 -hour  periods  (Group  VII)  at  all  horizons.  At  the  same  time 
at  the  station  located  6 miles  off  the  coast,  the  components  with 
this  period  did  not  exist.  However,  at  this  same  station  at  all 
depths  there  existed  components  with  periods  of  19.5  - 21.5  hours, 
which  were  not  revealed  at  the  offshore  station.  It  is  very  possible 
that  these  two  groups  of  components  were  generated  by  one  and  the 
same  cause,  i.e.,  from  the  diurnal  cycle  of  change. 

Apparently,  the  phenomenon  of  a diurnal  wave  is  transformed  and 
its  period  becomes  equal  to  19.5  - 21.5  hours  according  to  its 
proximity  to  the  coast.  In  our  opinion  it  seems  a probable  hypo- 
thesis that  causes  of  differences  in  and  the  absence  of  diurnal  com- 
ponents at  the  inshore  station  and  the  absence  of  components  with 
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periods  of  19.5  - 21.5  hours  at  the  offshore  station  are  explained  by 
the  disposition  of  the  stations  according  to  their  relationship  to  the 
nodes  of  standing  oscillations  with  analogous  periods. 

Fran  the  analysis  of  the  preceding  group  of  components  with  long 
periods,  there  is  likewise  clearly  seen  a distortion  (decreasing)  of 
periods  for  the  approach  of  waves  to  the  shore. 

As  far  as  the  series  of  observations  at  both  stations  are  concerned, 
the  internal  wave  periods  did  not  exceed  5 days;  the  maximum  possible 
periods  were  resolved  as  2-  and  3-day  components.  However,  it  is 
necessary  to  note  that  the  accuracy  of  these  determinations  was  very  low. 

All  enumerated  groups  of  oscillations  in  the  observed  inland  basin 
have  the  character  of  seiche  oscillations  and  a meteorological  origin  as 
the  basic  cause.  For  all  more  or  less  significant  atmospheric  perturba- 
tions (expressed  in  the  changing  plot  of  pressure  and  wind),  water  layers 
are  subjected  to  an  impulse,  and  oscillations  may  arise  in  them  with  the 
period  (cycle)  of  the  changing  pressure  and  wind. 

In  Fig.  2 there  are  presented:  (l)  a plot  of  pressure  and  wind 
(the  upper  curves)  at  the  coastal  station  at  the  time  the  multi-day 
stations  were  produced  and  for  the  12  preceding  days  (Fig.  2a);  (2)  a 
plot  of  pressure  at  the  offshore  station  (Fig.  2b);  and  (3)  curve  of 
sea  level  according  to  hourly  tide  gauge  readings  at  the  given  shore 
station  (Fig.  2c). 

As  clarified  by  the  analysis  of  surface  weather  charts  and,  in 
some  measure,  from  the  presented  patterns,  in  an  average  25-day  period 
over  the  region,  a deep  cyclone  (1006  m6\  ) is  found.  In  the  immediately 

succeeding  days  there  occurred  a rare  change  in  the  synoptic  situation. 
There  was  a cycle  span  of  high  pressure  at  this  place.  It  v/as  virtually 
destroyed  in  one  day.  From  the  27th  to  the  2nd,  a series  of  cycles 
traveled  across  this  region  at  intervals  on  the  order  of  2 days.  On 
the  3rd  there  was  produced  an  unusual  increase  of  pressure  for  records 
collected  in  the  region  of  the  station  at  the  crest  of  the  high  pressure, 
v/hich  to  the  end  of  the  3rd  formed  in  a nucleus.  In  the  following  days, 
synoptic  observations  determined  the  intensity  of  the  nucleus  and  the 
influence  of  the  trough  of  low  pressure,  located  to  the  northwest  of 
the  region.  By  observing  Fig.  2,  the  presence  in  the  plot  of  pressure 
and  wind  is  to  be  noted,*  also,  the  oscillation  of  sea  level  of  3 -day 
periodicity  is  a condition  due  to  the  strong  influence  of  the  nucleus 
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of  high  pressure  or  the  trough  of  low  pressure.  The  wind  in  the  region 
of  the  stations  and  at  the  nearby  shorepoint  was  very  weak  and  variable 
during  the  observation  period. 

It  is  quite  obvious  that  two-day  and  three-day  components  of 
internal  waves  in  the  series  taken  were  produced  from  an  exchange  of 
atmospheric  impulses.  Along  with  important  atmospheric  perturbations 
stimulating  long -period  oscillations,  in  all  stretches  of  the  pressure 
curves  viewed  there  were  small  perturbations  with  a variety  of  periods, 
which  nonetheless  may  have  initiated  seiche  oscillations  with  periods 
noted  in  the  tables. 

It  is  very  well  shown  that  such  agreement  may  serve  as  vhe  correla- 
tive tie  between  time  plots  of  atmospheric  pressure  and  density  at  the 
25 -meter  horizon  for  the  analyzed,  synchronized  stations.  Computation 
of  correlation  coefficients  were  made  by  the  formula: 


*(*)- 


N-t  V 
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rt  = l rtc=l  t-|»i 

~n=x  : 7*=?  \»  ir  n ; 
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where  .v  is  the  number  of  data  points,  -c  is  the  deity  time,  P is  the 
pressure,  and  a is  the  density  condition.  It  seems  that  for  Station 
II*  (inshore)  the  largest  correlation  coefficient  equaled  0.95  for  a lag 
time  of  nearly  29  hours.  At  the  same  time,  Station  II  (offshore)  had 
a maximum  correlation  coefficient  equal  to  0.34  without  delay  time. 


This  is  probably  a typographical  error  and  should 
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The  plots  of  pressure  at  both  stations  are  practically  one  and  the 
same,  and,  therefore,  for  the  attempt  of  constructing  the  correlative 
tie  between  density  oscillations  at  the  25-meter  horizon  at  both  stations, 
it  follows  to  anticipate  the  largest  correlation  coefficient  for  a lag 
time  on  the  order  of  29  hours.  Computations  of  such  correlative  ties 
show  that  the  greatest  correlation  coefficients  were  equal  to  0.71  for 
a delay  time  on  the  order  of  34  hours.  Moreover,  density  oscillations 
at  these  depths  in  the  region  of  the  inshore  station  were  delayed, 
relative  to  similar  oscillations  in  the  region  of  the  seaward  stations, 
by  34  hours.  It  was  possible  to  measure  the  speed  of  propagation  of 
the  observed  waves  by  knowing  their  delay  times.  Assuming  that  the 
wavefront  propagates  normal  to  the  coast,  its  speed  would  equal  36 
cn/sec;  if  the  front  propagates  at  an  angle  less  than  90°,  then  the 
speed  of  the  internal  wave  would  be  less  than  36  cm/ sec . 

For  the  isolation  of  periods  of  the  basic  components  of  the  waves, 
besides  the  formerly  stated  method,  we  used  the  Fourier  method  of  auto- 
correlation. That  method  isolates  the  components  in  order  of  their 
significance,  i.e.,  by  the  value  of  their  amplitudes  from  a series  of 
observations.  Formerly  / 1_/  we  used  the  contents  of  this  briefly 
stated  method. 

Periods  obtained  by  one  or  the  other  method  we re  found  to  be  in 
sufficiently  good  agreement.  V/e  will  examine  the  following  example. 

On  the  100 -meter  horizon  at  the  offshore  station  the  Fourier  analysis 
gave  the  value  of  the  period  of  the  predominant  component  as  26.9 
hours  (amplitude  6.3  meters),  and  the  second  most  significant  component 
had  a period  of  nearly  80  hours.  The  analogous  period  acquired  by  the 
second  method  were  equal  to  25.8  and  69  hours. 

For  the  50-meter  horizon,  t’ie  principal  period  acquired  was  equal 
to  25.7  hours  (amplitude  2.2  meters)  by  lie  first  method,  24.2  hours 
for  the  second  method. 

By  an  analogous  comparison  with  results  at  300  and  422  meters  of 
the  inshore  Station  I,  it  was  sho\/n:  (1)  the  Fourier  method  determined 
the  period  of  the  basic  components  at  the  390  -meter  horizon  — 21.1 
hours  (16.9  meter  amplitude) , and  at  the  422 -meter  horizon— 20.4  hours 
(12.8  meter  amplitude);  (2)  comparable  periods  acquired  by  the  applies.1  - 
tian  of  the  spectral  density  formula  were  20.3  hours  (300-meter  depth) 
and  21.1  hours  (422-meter  depth). 
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Therefore,  it  must  be  noted  that  there  are  good  comparative  re- 
sults acquired  by  both  methods. 

In  the  majority  of  cases  /"l,  5 _7,  application  of  the  Fourier 
method  gives  good,  reliable  results.  However,  when  the  predominant 
period  (according  to  amplitude)  of  the  next  largest  component  exceeds 
the  length  of  the  time  series,  or  when  the  time  series  is  less  than 
\ one -half  to  two  times  the  length  of  the  period,  computations 
by  the  Fourier  method  appear  difficult  for  the  determination  of  the 
periods  of  the  remaining  components . Such  difficulties  emerged  in 
our  computations  of  periods  for  the  100,  150,  and  200  meter  horizons 
at  Station  I (inshore)  by  these  methods.  By  the  use  of  the  computa- 
tional method  for  spectral  density,  these  difficulties  are  not 
encountered,  and  v/e  acquired  all  possible  periodic  components  existing 
at  the  time  of  observations.  However,  it  is  noted  that  the  long  periods 
acquired  were  quite  approximate  due  to  the  insufficient  extent  of  the 
stations . 

For  suitable  conditions  we  considered  it  expedient  to  use  either 
of  the  methods  for  the  acquisition  of  the  frequency  distribution  of 
internal  wave  components  in  similar  regions  of  the  sea. 

In  addition  to  spectra  of  internal  wave  periods  at  various  depths, 
periods  of  changes  in  atmospheric  pressure  and  sea  level  were  also 
determined.  The  data  obtained  are  reduced  in  Table  3.  As  a result  of 
the  comparisons  in  the  three  tables,  there  is  an  obviously  good  rela- 
tionship between  periods  of  oscillation  on  all  the  previously  mentioned 
horizons  and  periods  of  sea  level  and  pressure  changes.  The  results 
of  the  comparisons  once  more  confirm  the  opinion  that  in  deep  inland 
seas  atmospheric  impulses  are  shown  to  influence  water  layers  very 
significantly . 

On  the  basis  of  all  the  foregoing  one  nay  make  the  following 
conclusions : 

1.  Each  of  the  observed  curves  of  density  isopleths  (for  the 
offshore  as  well  as  the  inshore  station)  appears  in  the  nature  of  a 
complex  oscillation,  consisting  of  a sum  of  vertical  oscillations  of 
various  periods  (as  in  Tables  1 and  2) . 
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oscillations  appeared  to  be  dominated  by  components  with  periods  of 
21-24  hours  and  periods  on  the  order  of  3 days. 

3.  The  principal  factor  in  the  phenomena  of  internal  waves 

in  deep  inland  seas  appears  to  be  the  variability  of  atmospheric  condi- 
tions, e.g.,  pressure,  wind  regime,  etc.  In  the  cases  observed,  because 
of  the  effect  of  the  particular  positions  of  the  stations  in  their 
proximity  to  the  center  of  the  barometric  phenomena,  observed  winds 
were  very  weak  and  variable  in  direction.  Accordingly,  in  the  given 
case,  the  principal  v/ave  excitation  factor  appeared  as  the  variability 
of  atmospheric  pressure. 

4.  The  semidiurnal  tidal  component  exists  at  all  horizons  at 
both  stations  and,  in  spite  of  this,  it  does  not  appear  to  be  dominant 
but  has  a large  significance  in  the  mixing  of  water  at  great  depths. 

5.  For  the  objective  estimate  of  phenomena,  it  is  necessary 
to  limit  the  isolation  to  one  or  two  components  of  the  waves.  The 
spectral  analysis  of  internal  waves  makes  it  possible  to  determine  the 
degree  of  influence  of  one  or  another  of  the  forces  exciting  internal 
waves  with  corresponding-  parameters.  Knowing  the  state  of  these  forces 
and  the  laws  governing  them,  it  may  further  be  possible  to  construct 
laws  about  the  propagation  of  internal  waves.  Finally,  as  noted  earlier, 
this  knowledge  has  great  practical  value. 

The  study  of  the  oscillations  of  internal  seiches  in  deep  inland 
seas  is  of  interest  because  water  movement  in  the  vertical  and  hori- 
zontal directions  at  great  depths  is  accompanied  and  caused  by  them. 

For  the  acquisition  of  more  complete  and  reliable  data  for  the 
quantitative  expression  of  the  components  of  internal  waves,  it  is 
necessary  to  occupy  multi -day  stations  for  a sufficient  duration  not 
only  inshore  but  in  deep  water  parts  of  the  sea  and  with  the  partici- 
pation of  several  ships  working  in  synchronization  according  to  a 
special  program. 

In  conclusion,  the  authors  express  deep  gratitude  to  A.  G. 
Kolesnikov  for  valuable  discussions  and  help  in  the  work,  and  also 
A.  F.  Sokonnov,  who  conscientiously  executed  a valuable  part  in  Uie 
laborious  computations. 
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Fig.  1 - Distribution  of  spectral  density  on  the  300-meter  horizon  at  the  station  closest  to  «»w 
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Fig.  2 - Changes  with  time  in  meteorological  elements  and  sea  level  according  to  date  from  coastal 
and  open-sea  stations 

a.  Change  in  pressure  and  wind  at  coast  station 

b.  Change  in  pressure  at  sea  station 

c.  Change  In  sea  level  at  shore  station 
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